Search for: All records

Abstract Regularization of seismic inversions has a strong imprint on tomographic images. We analyze recorded and spectral‐element S, Sdiff, and SS waveforms to evaluate the benefit of body‐wave amplitudes in global tomography. L‐curve analysis for S40RTS models with recorded and synthetic waveforms show that SS‐S traveltimes and SS/S amplitude ratios have minima within the same damping parameter range. SS/S ratios for S40RTS and model GLAD‐M25 show the trade‐off between scale‐length and strength of lowermost‐mantle heterogeneities. The recorded SS/Sdiff ratios are lower than predicted by 3D mantle models which may be explained by a decrease in the mean shear velocity by at the lowermost 200 km of the mantle. Our results suggest that SS/S amplitude measurements made with 3D waveforms can be used to constrain damping in linearized inversions, and amplitudes are essential for studying the size of heterogeneities. 

Stress drop, a crucial source parameter in earthquake studies, significantly influences ground motion prediction and seismic hazard assessment. Despite several existing methods to estimate stress drops, the resulting stress drop estimates often exhibit a wide variation of up to 3-4 orders of magnitude. In this study, we address the robustness of stress drop estimation by introducing a point-wise spectral ratio stacking approach based on empirical Green’s functions (eGfs). Conventional trace-wise stacking can lead to data exclusion due to high signal-to-noise ratio requirements across a wide range of frequency. By adopting point-wise stacking, we maximize the utilization of useful recording information, leading to more accurate stress drop estimates. We applied the point-wise spectral ratio stacking method to a comprehensive dataset comprising global earthquakes from 1990 to 2020 with magnitude larger than Mw5.5 and depth shallower than 50 km. We first verified the moment magnitudes of earthquakes estimated from the resulting seismic moment ratios. We found that the moment magnitude of master events best consistent with catalog magnitudes when the magnitude difference between master and their eGfs differs by about 0.5. Our analysis indicates that stress drop of shallow earthquakes exhibits no depth dependence, while showing a slight increase with magnitude. The results obtained through our optimized stacking process shed new light on stress drop estimate of shallow earthquakes and have the potential to enhance the understanding of earthquake mechanics. 

ABSTRACT Seismic tomography is the most abundant source of information about the internal structure of the Earth at scales ranging from a few meters to thousands of kilometers. It constrains the properties of active volcanoes, earthquake fault zones, deep reservoirs and storage sites, glaciers and ice sheets, or the entire globe. It contributes to outstanding societal problems related to natural hazards, resource exploration, underground storage, and many more. The recent advances in seismic tomography are being translated to nondestructive testing, medical ultrasound, and helioseismology. Nearly 50 yr after its first successful applications, this article offers a snapshot of modern seismic tomography. Focused on major challenges and particularly promising research directions, it is intended to guide both Earth science professionals and early-career scientists. The individual contributions by the coauthors provide diverse perspectives on topics that may at first seem disconnected but are closely tied together by a few coherent threads: multiparameter inversion for properties related to dynamic processes, data quality, and geographic coverage, uncertainty quantification that is useful for geologic interpretation, new formulations of tomographic inverse problems that address concrete geologic questions more directly, and the presentation and quantitative comparison of tomographic models. It remains to be seen which of these problems will be considered solved, solved to some extent, or practically unsolvable over the next decade. 

ABSTRACT Although the Brune source model describes earthquake moment release as a single pulse, it is widely used in studies of complex earthquakes with multiple episodes of high moment release (i.e., multiple subevents). In this study, we investigate how corner frequency estimates of earthquakes with multiple subevents are biased if they are based on the Brune source model. By assuming complex sources as a sum of multiple Brune sources, we analyze 1640 source time functions of Mw 5.5–8.0 earthquakes in the seismic source characteristic retrieved from deconvolving teleseismic body waves catalog to estimate the corner frequencies, onset times, and seismic moments of subevents. We identify more subevents for strike-slip earthquakes than dip-slip earthquakes, and the number of resolvable subevents increases with magnitude. We find that earthquake corner frequency correlates best with the corner frequency of the subevent with the highest moment release (i.e., the largest subsevent). This suggests that, when the Brune model is used, the estimated corner frequency and, therefore, the stress drop of a complex earthquake is determined primarily by the largest subevent rather than the total rupture area. Our results imply that, in addition to the simplified assumption of a radial rupture area with a constant rupture velocity, the stress variation of asperities, rather than the average stress change of the whole fault, contributes to the large variance of stress-drop estimates. 

Abstract The upper boundary of the mantle transition zone, known as the “410-km discontinuity”, is attributed to the phase transformation of the mineral olivine (α) to wadsleyite (β olivine). Here we present observations of triplicated P-waves from dense seismic arrays that constrain the structure of the subducting Pacific slab near the 410-km discontinuity beneath the northern Sea of Japan. Our analysis of P-wave travel times and waveforms at periods as short as 2 s indicates the presence of an ultra-low-velocity layer within the cold slab, with a P-wave velocity that is at least ≈20% lower than in the ambient mantle and an apparent thickness of ≈20 km along the wave path. This ultra-low-velocity layer could contain unstable material (e.g., poirierite) with reduced grain size where diffusionless transformations are favored. 

SUMMARY Long-period (T > 10 s) shear wave reflections between the surface and reflecting boundaries below seismic stations are useful for studying phase transitions in the mantle transition zone (MTZ) but shear-velocity heterogeneity and finite-frequency effects complicate the interpretation of waveform stacks. We follow up on a recent study by Shearer & Buehler (hereafter SB19) of the top-side shear wave reflection Ssds as a probe for mapping the depths of the 410-km and 660-km discontinuities beneath the USArray. Like SB19, we observe that the recorded Ss410s-S and Ss660s-S traveltime differences are longer at stations in the western United States than in the central-eastern United States. The 410-km and 660-km discontinuities are about 40–50 km deeper beneath the western United States than the central-eastern United States if Ss410s-S and Ss660s-S traveltime differences are transformed to depth using a common-reflection point (CRP) mapping approach based on a 1-D seismic model (PREM in our case). However, the east-to-west deepening of the MTZ disappears in the CRP image if we account for 3-D shear wave velocity variations in the mantle according to global tomography. In addition, from spectral-element method synthetics, we find that ray theory overpredicts the traveltime delays of the reverberations. Undulations of the 410-km and 660-km discontinuities are underestimated when their wavelengths are smaller than the Fresnel zones of the wave reverberations in the MTZ. Therefore, modelling of layering in the upper mantle must be based on 3-D reference structures and accurate calculations of reverberation traveltimes. 

SUMMARY A number of seismological studies have indicated that the ratio R of S-wave and P-wave velocity perturbations increases to 3–4 in the lower mantle with the highest values in the large low-velocity provinces (LLVPs) beneath Africa and the central Pacific. Traveltime constraints on R are based primarily on ray-theoretical modelling of delay times of P waves (ΔTP) and S waves (ΔTS), even for measurements derived from long-period waveforms and core-diffracted waves for which ray theory (RT) is deemed inaccurate. Along with a published set of traveltime delays, we compare predicted values of ΔTP, ΔTS, and the ΔTS/ΔTP ratio for RT and finite-frequency (FF) theory to determine the resolvability of R in the lower mantle. We determine the FF predictions of ΔTP and ΔTS using cross-correlation methods applied to spectral-element method waveforms, analogous to the analysis of recorded waveforms, and by integration using FF sensitivity kernels. Our calculations indicate that RT and FF predict a similar variation of the ΔTS/ΔTP ratio when R increases linearly with depth in the mantle. However, variations of R in relatively thin layers (< 400 km) are poorly resolved using long-period data (T > 20 s). This is because FF predicts that ΔTP and ΔTS vary smoothly with epicentral distance even when vertical P-wave and S-wave gradients change abruptly. Our waveform simulations also show that the estimate of R for the Pacific LLVP is strongly affected by velocity structure shallower in the mantle. If R increases with depth in the mantle, which appears to be a robust inference, the acceleration of P waves in the lithosphere beneath eastern North America and the high-velocity Farallon anomaly negates the P-wave deceleration in the LLVP. This results in a ΔTP of about 0, whereas ΔTS is positive. Consequently, the recorded high ΔTS/ΔTP for events in the southwest Pacific and stations in North America may be misinterpreted as an anomalously high R for the Pacific LLVP. 

Seismology provides important constraints on the structure and dynamics of the deep mantle. Computational and methodological advances in the past two decades improved tomographic imaging of the mantle and revealed the fine-scale structure of plumes ascending from the core-mantle boundary region and slabs of oceanic lithosphere sinking into the lower mantle. We discuss the modeling aspects of global tomography including theoretical approximations, data selection, and model fidelity and resolution. Using spectral, principal component, and cluster analyses, we highlight the robust patterns of seismic heterogeneity, which inform us of flow in the mantle, the history of plate motions, and potential compositionally distinct reservoirs. In closing, we emphasize that data mining of vast collections of seismic waveforms and new data from distributed acoustic sensing, autonomous hydrophones, ocean-bottom seismometers, and correlation-based techniques will boost the development of the next generation of global models of density, seismic velocity, and attenuation. ▪  Seismic tomography reveals the 100-km to 1,000-km scale variation of seismic velocity heterogeneity in the mantle. ▪  Tomographic images are the most important geophysical constraints on mantle circulation and evolution.